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INTRODUCTION 


A nupt>er  of  techniques  have  been  developed  for  modeling  the  propagation 
of  ultrasonic  waves  in  layered  media, and  a number  of  authors  have  also 
discussed  the  case  for  which  such  media  are  infinitely  periodic  in  spaceU-^.5. 
However,  none  of  these  authors  has  fully  treated  the  case  of  reflection  and 
transmission  through  finite  periodic  media.  Brillouin'^  has  indicated  these 
cases  might  be  treated  whenever  boundary  conditions  are  such  that  no  waves  are 
reflected  from  the  e.xtremity  of  the  periodic  medium,  but  has  not  worked  out 
this  problem  in  detail. 


In  this  report  a general  technique  for  treating  arbitrary  finite  periodic 
media  for  all  boundary  conditions  is  derived  and  the  interesting  case  in  which 
the  ambient  media  at  the  two  boundaries  have  identical  acoustic  properties  is 
discussed  in  detail. 


This  generalization  is  accomplished  by  taking  the  iterative  equations 
derived  in  reference^  and  reformulating  them  in  such  a manner  as  to  represent 
each  iteration  as  matri.x  multiplications.  A matrix  diagonalization  technique 
is  then  used  to  produce  a closed  form  expression  for  an  arbitrary  number  of 
iterations  "q"  corresponding  to  a laminate  of  "q"  layers. 

The  resulting  formulae  rigorously  confirm  many  of  the  empirical  obser- 
vations of  reference^  relating  to  the  spacings  of  thickness  resonances  and 
the  positions  of  forbidden  frequency  bands. 


THEORY 

A.  Wave  Propagation  in  Layered  Media 

Techniques  have  been  reported^  for  calculating  the  reflection  and  trans- 
mission coefficient  of  layered  media  when  the  transmission  and  reflection 
coefficients  at  each  of  the  interfaces  are  known. 

These  equations  may  be  applied  to  an  array  of  n laminae  of  the  type  shown 
in  figure  1.  The  medium  to  the  right  of  the  first  interface  is  assumed  to  be 
infinite  in  extent  as  is  the  medium  to  the  left  of  the  n^^i  interface.  In 
general  a given  medium  i has  a thickness  d^,  a velocity  Cj^  and  has  reflection 
and  transmission  coefficients  ri  and  ti  associated  with  waves  incident  from 
the  left  upon  its  left  interface  and  has  frequency  dependent  reflection  and 
transmission  coefficients  Rj^  and  Tj^  associated  with  the  entire  stack  of 
laminae  to  the  right  of  that  interface. 

Applying  the  equations  of  reference  1 to  the  laminate  of  figure  1 gives 
the  following  results: 


* 

S'-  BUB^not 
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Rp  = ro  Tq  = tp 

(1) 

p , . _ ri  + ro  exp(-2ioodi/Ci) 

^ 1 + rj  ro  exp(-2iwd}/C^) 

(2) 

. . tl  ^0  exp(-iudi/Ci) 

* 1 + rj  ro  exp(-2iwdi/Ci) 

(3) 

and  in  general 

p C ■)  _ Kn-i(w)  exp(-2iwdn/Cn) 

" “ 1 ♦ ^n  Rn-i(w)  exp(-2icodn/Cn) 

(4) 

T (co)  = exp(-ia)dn/Cn) 

" 1 l<n-lC“)  exp(;-2ia)dn/Cn)  . 

In  order  to  make  these  equations  more  tractable  it  is  desirable  to  transform 
them  from  an  iterative  into  an  explicit  form;  therefore,  the  following 
definitions  are  made: 


^n  = >"n  ♦ Kn-1  exp(2i()>n)  (6a) 
i>n  ■ "wdn/Cn  (6b) 
^’n  = 1 + Rn-i  exp(2i()>n)  (6c) 
from  which  it  follows  that 

but 


Rn  = 


rn  bn_i  ^ exp(2i»n)  a„_i 

b^-l  ♦ rn  exp(2i(|>n)  an_i 


(8) 


Clearly  the  two  quantities  an  and  bn  may  then  be  written  as  a vector  and 
equation  8 becomes: 


\Ki 


exp(2i4in)  tn 
rncxp(2i4>n)  1 bn.i 


(9) 


B.  Periodic  Layered  Media 
1.  Reflectance  liquations 

Since  only  the  ratio  an/bn  is  of  interest  a more  symmetrical  form  may  be 

used ; 
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I exp(i4>n) 
rncxp(i<j>n) 


rj^exp(-i(j)n) 

exp(-i()!„) 


Si-1 

3n-l 


no) 


which  immediately  implies, 
<^xp(i4>m) 
r„exp(i<lin,) 

m=  1 , n 


1 

7T  1 

1 Bn  I ’ 

\ 

\ 

r„,exp( \ 
exp(-i(j>ni) 


/J 


ro 

1 


(11) 


ITiis  form  is  somewhat  suggestive  of  the  result  of  Brekhovskikh^;  however, 
it  has  the  advantage  that  it  contains  fewer  implicit  operations. 

For  periodic  structures  expression  (11)  takes  the  form. 


\i 

Bn 


7T 

m=l , ^ 


exp(i()>m)  rn,exp(-i<frm)  ^ 

rmexp(i<l>m)  exp(-i4>n,)^ 


(12) 


1 ^91^ 

1 ’’0  ^ 

1 ^ 

\^9il 

1 = Mp'’ 

; where  q = 

(12a) 

^ ^ i 

n is  the  total  number  of  layers  in  the  structure,  is  the  number  of 
layers  in  one  period  of  the  structure^  IL  is  the  number  of  periods  in  the 

I 

laminate,  and  Mp  is  the  product  matrix  in  expression  12. 

Since  the  product  matrix  in  (12)  will  not  in  general  be  hermitean,  an 
expression  for  a 2x2  matrix  raised  to  a power  is  required.  For  any  matrix 
M it  can  readily  be  shown  that  for 


Mp  = 


f 

h 


Mp  • Mp  = (e+h)  M + [2+fg-eh]  1 
where  11  is  the  unit  matrix 


(13) 

(14) 


Since  multiplying  M by  a constant  factor  will  not  affect  the  value  of 
Mp  may  be  used. 

det (Mp) 


the  matrice  M 


M • M = 2aM*  1 
e+h 


(15) 


where  o = 


2(fg-oh) 


1 
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It  is  clear  from  (15)  that  the  product  fg-eh  must  be  of  the  form 

(16) 

This  may  be  evaluated  for  a specific  case  by  repeated  application  of 
expression  (15)  or  by  using  the  matrix  relation 


(17) 


The  2x2  matrix  in  (17)  is  s>Tnmetric  and  has  eigenvalues  X = a ± a^+l 
and  eigenvectors  (1,  a + a^+l)  and  (1,  -a  - a^+1).  Then  for  a real 

this  matrix  may  be  diagonalized  as  follows: 


/ ““N  \ / 

N-1 

f “0  ^ 

1 2- 

1 \ N-1 

[ 1 1 1 

" i 

i j 

i ^ 

0 1 

rl 

/2a  0 

N -(a2.1)N 

1 1 oj 

(a  -A/a^+l) 

(a-V^ a^  + 1)^*^  - (a+A/a^  + l ) ^ ^ + - (a+\/a^  + 1)^ 


(a-Va^  + l)^  - (a+Va^  + 1)  ^ (a ^/a^  + l )^” ^ - («A/a^  + l) 


iN  (a"  + l) 


M'"  - 


N 


(a  - a/^)  ^ ( a + V^)  1 


(a  - "Va^* 1 ) 

2.  Travelling  Wave  Regime 


M+n  j (a- Va^  + l)’^-(ct+Va^l)*''| 

(18)- 


As  an  alternative  to  applying  expression  (18)  we  may  diagonalize  the 
matrix  in  12a  directly.  The  form  of  this  matrix  is  easily  shown  to  be, 

★ 


M^l  M21 
M21  Mjj 


(19) 


In  the  most  general  case  such  a matrix  will  have  two  distinct  eigenvalues 
>♦  given  by 


X±  = 1/2 


Mil  * ^11  l.~y^(Mi  ^+Mj  1)  ^ - 4 [ I Mj  j|  I M21I  ‘■] 


for  (Mj,  ♦ m|i)2  < 4 [IMiil^  + IM2ll^]. 

which  is  the  range  of  periodic  solutions,  X±  take  the  form 
Ai  = P 

Mil  * Mil 


where  cos  0 = 


2p 

2 ^ I u_ - I 2 1 


and  p7  = [ |Mi  il  ♦ I M21 1 

and  from  (12)  p^  = = 1,11  (1-r^^) 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 

(20f) 
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Corresponding  to  these  two  eigenvalues  two  non-orthogonal  eigenvectors  may  be 
constructed  of  the  form 


V* 


(21) 


Suppressing  the  I subscript  12a  may  be  written. 


ype'^^'^V+  ♦ vpe'^p'^V- 


where  u and  v are  the  components  of  the  vector 
independent  basis  vectors  V+  and  V-. 


with  respect  to  the 


(22) 


Again  as  in  (10),  ignoring  extraneous  proportionality  factors  Rq  may  be 
expressed  in  the  form 


where 

A,n 

'-A'’+M2j  -A‘’-M22  \ 

Mj  2 ■ A M2^ 

, ^‘’♦(Mii-AJ  A^_(Mii-Aj/ 

(23) 


In  order  to  illustrate  the  nature  of  one  such  solution  let  us  consider 
the  important  case  suggested  by  Brillouin  in  which  tq  = 0,  This  may  be 
realized  by  placing  an  infinite  medium  on  the  right  side  of  the  array  which 
is  identical  to  the  last  medium  on  the  right  side  of  the  array  (see  figure  1) . 


Solving  for  the  coordinate  in  the  eigenvector  space  gives  the  simultaneous 
equation 


-p  ^2 2 "V  ^21  ~ ^ 

u(Mh-A^)  + v(Mji-AJ  = 1. 


Which  imply 

V = -p 

p = [A.-Aj"*  = (2pisin\|)]  . 

In  fact  the  value  of  p is  of  no  consequence  since  it  cancels  in  the  final 
expression  and  we  have 
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R 


q 


-j\\  - X^)  M21 


Mt2sinq\j^ 

-Mijsinqi|)  + p^^'^(qH)ii^ 


R 


q 




pSin^cotqiJj  - (Mj ^ -p‘“°^i|() 


(24a) 


(24b) 


Since  Mjj  and  Mj2  are  sums  of  products  of  bounded  periodic  functions  R will 
will  have  zeroes  wherever  cot  ij;  -*•  ® 


i);  = — for  n an  integer  / 0 
q 


(25) 


Then  assuming  that  i(<  is  a monitor  function  of  frequency  Rq  will  have  q -1 
minima  for  which  Rp=0  between  \J;=u)=0  and  ^=‘n . For  cases  in  which  — is  an 
integer  (24b)  has  an  indeterminant  form  in  the  denominator  and  each 
case  must  be  considered  carefully. 


In  section  (5)  it  is  demonstrated  that  this  result  may  also  be  obtained 
for  a more  general  set  of  boundary  conditions  in  which  the  ambient  media 
bordering  the  array  from  the  right  and  left  are  identical. 


3.  Non-Travelling  Wave  Solutions 

For  (Mil  * Mii*)2  > 4 ()Mu|2  + |M2i|2) 
again  set  p2  = [|Mn|2  + |M2i|2] 


and  X^  = 
where  y+ 


Pp 


yt 


I 


= loge 


hl‘*'^ll  ^'^(Mii+Mii  )^-4p 


It  is  a straightforward  substitution  to  show  that  e^'*’  = (e^*)'^ 
Then  letting  y = y*  the  expression  for  the  eigenvalues  becomes, 


X+  = p^/2,±Y 


(26) 


Together  with  (24a)  this  gives 


- (sinh(qY))  M21 

Rq  

Mil (sinh(qY))  - (sinh(q+i)Y) 


Kq  = ^ (27) 

pcoshY  + psinhy  cothY  - Mu 
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The  hyperbolic  functions  are  not  periodic  and  over  a range  of  monoto- 
nically  increasing  y they  should  exhibit  at  most  one  maximum  or  minimum. 

4.  Phase  Velocity 

The  phase  velocity  of  sound  for  a finite  plate  in  a uniform  ambient 
medium  has  been  shown  to  he. 


c: 


2df 

n 


(28) 


where  f is  the  frequency  of  the  n^^  thickness  resonance  from  a plate 
laminate  of  thickness  d. 


It  will  now  he  ’^hown  that  (28)  when  used  as  a definition  for  phase 
velocity  allows  a generalization  of  this  concept  to  periodic  layered  media. 

It  is  by  no  means  intuitively  obvious  that  phase  velocity  in  a laminate  is 
a material  property  uniquely  defined  by  (28),  since  the  positions  of  thick- 
ness resonances  could  be  affected  by  boundary  conditions  in  a thick, .ess 
dependent  manner.  This  in  fact  appears  to  be  the  case  whenever,  there 
exists  at  the  far  interface,  a non  zero  reflection  amplitude  off  a medium 
which  is  different  from  the  incident  medium.  Clearly  this  limits  the 
measurement  situations  in  which  the  definition  (28)  can  be  applied.  The 
definition  is  further  limited  in  that  it  can  only  be  applied  to  the  deter- 
mination of  velocity  at  particular  frequency  intervals.  However  this  should 
not  limit  the  validity  of  the  definition  as  long  as  it  is  not  inconsistent 
for  laminate-  of  different  thicknesses.  It  will  now  be  demonstrated  that 
thi-  IS  in  fact  true  whenever  the  relations  (25)  and  (20d)  are  satisfied. 

from  equation  (24b)  such  frequency  resonances  occur  in  periodic  lam- 
inates whenever 

nw 

^ = T (25) 


where  is  implicitly  defined  by  the  relationship 


cos^  = ^11  * 


2p 


(20d) 


Since  (25)  and  (28)  are  both  functions  of  1 n | and  therefore  not  dependent 
on  d and  since  iji  is  also  independent  of  d then, 


C = 


2v  Df 

*(f) 


uiD 

^Kfy 


(29) 


Ibis  yields  the  general  result  that  the  phase  velocity  of  sound  in  any 
periodic  medium  as  defined  by  (29)  is  independent  of  the  thickness  of  that 
medium  whenever  (25)  and  (20d)  are  satisfied.  The  next  section  will  explore 
some  general  conditions  for  which  these  relations  hold. 
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5.  The  Effects  of  an  Ambient  Medium 


In  sections  (2)  and  (4)  above  an  expression  was  derived  for  reflectance 
and  phase  velocity  for  any  periodic  media  for  which  the  acoustic  properties 
of  the  outer  boundary  layers  match  those  of  the  ambient  medium.  In  this 
section  we  will  generalize  those  results  to  periodic  media  for  which  the 
acoustic  properties  of  the  outer  boundary  layers  are  identical  to  each  other 
but  not  necessarily  to  those  of  the  ambient  medium.  This  case  is  approxi- 
mated by  any  composite  laminate  emersed  in  a water  bath  for  ultrasonic 
interrogat ion . 


In  general  the  matrix  representing  such  a situation  can  be  written  in 
the  form 


M. 


= MsM3^Mp)"Ms  ^ 


0 


/ 


(30) 


where  Mp  is  the  product  matrix  for  one  laminate  period 
q is  the  number  of  laminate  periods 

Mq  is  the  matrix  for  the  outer  layer  at  the  laminate  interior 

Mg  is  the  matrix  for  the  outer  layer  at  the  incident  interface 

Mg  is  the  matrix  for  an  arbitrary  layer  of  incident  medium. 


(30)  can  readily  be  transformed  into  the  form, 
/ 


Mp 

i'l'l 
e i4'l 


.1' 


1 )=  MgMo'^Mg  (Mg-lMpMg)^ 


-i^l 


■i(J>l 


.i'J'l, 


’ll  '’21 


\ ^21 


M 


11  / 


' 1 y' 


Y 1 


/J 


where  <ij  is  the  phase  shift  associated  with  the  outer  layer  and  y is 
the  reflection  amplitude  coefficient  for  a wave  reflected  from  the  outer 
medium. 


(1-Y^)‘^ 


Mii'Y^Mjj  y(Mi]-Mjj)  +M2i*-y^M2i 

y(Mj j*-Mjj)+M2i-Y^M2j  M^^*-Y^Mjj 


0 

1 


(31) 


which  from  20d  yields. 


COSli) 


M] j+Mjj 


2p 


(32) 


This  implies  that  the  phase  velocity  as  defined  in  section  (4) 
to  a large  class  of  periodic  layered  media. 


applies 
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C.  Two  Component  Media  fA  Composite  Approximation) 

1.  Reflectance  for  Travelling  Waves 

Consider  the  periodic  medium  of  figure  2.  Applying  expressions  (12)  and 
(12a)  gives  the  expression. 


/ ^ n \ 

e 0 

' e^'^’A 

1 0 e ‘^B/ 

-i<^A 

\ - r e 

' i (4>A+2(Jib)  2 

e -re 

i(24'B'‘t’A) 

, -i<liA  *'^A, 
r ( e e ) 

-(2<t>B'"<'A)_^2^i('('A-2<t>B) 


from  expression  20d 

cos(2()>g+(|)^)  - r2cos(2(J>g-()i^) 

cosil)  = 

(l-r2) 

From  24b  for  a laminate  of  q periods, 


(33) 


(34) 


irs  inij)^ 

psin<J^cotqi(/  - (g^'^A+2i(}iB_j.2  -pcostii) 


(35) 


2.  Velocity  of  Sound-Low  Frequency  Limit. 

The  use  of  expression  (29)  for  velocity  of  sound  requires  the  inversion 
of  (34)  which  produces  a rather  unwieldy  result;  however  in  certain  limits 
approximations  can  be  made  which  produce  tractable  expressions.  For  (|iy^+24>g 
small  (34)  becomes 


cosi;/  = 1 ♦ 


1/2  (A'l>g^^4((iy^4'g*4i/v^)-t^/2('*<t'B^‘‘*'*’A‘^B*‘*’B^^ 


(36) 


1 -r'^ 


1 ♦ ^ 


(2(>b<>a)‘ 


lOr'^cJi^^ip 
' l-r^" 


1/2 


(37) 


‘*A  ♦ ‘^B 


2 li  ^ ^A<^B 

Cb^Ca]  1-1-2  CaCb_ 


(38) 
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FIGURE  2.  Two  Component  Layered  Media 
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5.  \clocity  of  Sound  - Lirn  i t of  Thin  Boundary  Layers  and  Low  Frec^uenc'/ , 


Again  from  (34)  , 


costp 


cos  (2?g+ii^)  - r-  cosf2>g-}^) 
(1-r-) 


For  25B<<ii;^ 


COS'I)  = 


cos4i.^+2(sini)^)  ig 


(r--I) 

(r-+i) 


(391 


for  j)  small 


cosii 


cos (5^- 2 1 


1 r-  + l 


From  (30)  , 


ujD  _ u)  (d^-t-Zd^) 

jid^  2r-  2dBuj 

~c2  *(l-r-)  Cb 


C 


‘iA’ 


i J 


2dBI 


A I dA/Cg  l-r^\  dA/  C. 


dA*2dB 

^ ^ d^ 

Ca  Cb  ^ 


C40) 


The  third  term  in  the  denominator  of  l40)  describes  the  effect  of  inter- 
nal reflections  retarding  the  phase  of  the  wave  train  passing  through  the 
laminate.  The  expression  (40)  will  probably  be  useful  for  relating  static 
and  dynamic  (i.e.,  ultrasonic)  measurements  of  modulus  at  low  frequencies. 

When  r=0  the  internal  reflections  disappear  and  (40''  reduces  to  the  "a  priori" 
result 


C 


o 


dA+2dB 


'Ia  -^b 


iC 


(41) 


which  is  the  usual  expression  for  average  velocity.  It  is  of  interest  that 
both  (40)  and  (41)  are  dependent  only  on  the  ratio  of  dA/^dg  (i.e.,  the  percent 
composition)  and  not  on  specific  thickness.  For  this  reason  (41)  may  be  a 
reasonable  approximation  for  the  low  frequency  velocity  of  sound  in  materials 
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with  randomly  spaced  layers  such  as  directionally  solidified  eutectics. 

4.  Discontinuity  in  Phase  Velocity  for  a Two  Component  Laminate  with  Thin 
Boundary  Layers 


Prom  (341 


COSlJ/ 


cos(2<t>B+4>A)-r^  cos (2<}iB-4iy^) 


(l-r2) 

let  itiy^>>24iB  and 
Then  for  4'a~'^  we  have, 

, , 2 

COSlJ^  = -1+  


cos4> 


l-r2 

4’a^  2r2-t-l  <Ji 


A(}>B 


2 r2-l 

In  order  to  find  the  region  of  phase  velocity  discontinuity  we  set 

7 2r2+l 

*A<fB  * = 0 

r^-1  2 


(42) 


which  produces  the  condition 


= -2^  ^B 

r2-l 


1± 


r2.1 


1-r" 


-1 


(43) 


The  approximate  frequency  and  width  of  the  discontinuity  may  now  compute 
from  (43) 


r2+l 

r2-l 


= -2 


U = 


For  r=0  this  reduces  to  fg= 


2 

-1 

-V  |r2-lj 

2dA  4dB 

r2+l 

r^l  -1 

^A  ^B 

1 -r2 

-1 


(44) 


2tlA  ^ ^ 


1 -1 


which  is  the  "a  priori"  first 


approximation  to  the  thickness  resonance  frequency  for  the  two  layers  A and  B 
together  and  no  forbidden  gap  exists. 
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^ • 4 ’ 1 I 

j physically  allowable  values  of  r,  is  is  clear  that  the 

onset  of  the  first  forbidden  band  is  below  fg  and  since 


the  end  of  the  forbidden  band  will  always  be  below  the  thickenss  resonance  for 
the  acoustically  thickest  layer 


fOA 


(45) 


for  computational  purposes  when  dealing  with  small  discontinuities  it  is 
most  convenient  to  use  the  relationship, 


A 


1 

l+r^ 

(l+r2)\^  1 

- 8 — 
<^B 

l-r2 

‘J 

1/2 


which  is  derived  directly  from 


If  iUr) 


l*r^ 

1 -r2 

(l-r2)2 


-1 


(44). 
1 1/2 


(46)  becomes, 


f2 


8dc 


n(r) 


(46) 


(47) 


Since  fl(r)  varies  slowly  with  r (See  figure  3)  even  a very  rough  estimate 
of  r will  allow  a fairly  accurate  determination  of  the  bond  line  thickness 
(2dg).  Conversely  if  the  bond  line  thickness  is  known  with  reasonable 
accuracy,  a very  accurate  determination  of  the  value  of  r can  be  made. 


NUMERICAL  RESULTS 


Water/Class  Arrays 

Equation  (25)  is  an  important  result  since  it  allows  the  prediction  of 
ultrasonic  reflection  minima  for  any  periodic  medium  subject  to  the  conditions 
established  in  section  B-5  above.  This  prediction  requires  only  that  the  func- 
tion if/  be  tabulated  as  a function  of  ui. 

As  a numerical  check  on  the  validity  of  Equation  (25)  and  the  results  of 
section  B-5  the  reflectance  spectrum  for  an  array  of  water  and  glass  laminae  in 
in  an  ambient  medium  of  ethyl  alcohol  was  computed  using  the  formulae  (1) 
through  (5).  The  total  number  of  media  involved  were  11  and  these  constituted 
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a four  period  array  (i.e.,  q =4)  in  a non  matching  ambient  medium. 

Table  I compares  these  calculations  with  ones  obtained  from  equation  '25). 
Agreement  is  good  between  the  two  sets  of  calculations  e.xcept  for  the  case  n = 4 
which  IS  an  indeterminant  form  for  eq  (24b)  for  which  relation  (25)  is  no 
longer  valid.  In  principal  both  calculations  are  e.xact  and  should  give  iden- 
tical results;  however  because  of  the  large  number  of  iterations  (11)  required 
for  formula  (1)  through  i)5)  the  propagation  of  error  involved  with  the  trigo- 
nometric algorithms  probably  becomes  significant. 

Graphite/ Epo.xy  Laminates 

In  section  C-1  a general  relationship  is  derived  between  the  phase  and 
the  frequency  for  waves  reflected  from  a two  component  laminate  and  through 
the  use  of  e.xpression  (29)  this  may  be  converted  into  a dispersion  relation 
C=C(o.j.  Figure  4 is  a plot  of  such  a disoersion  relation  for  a hvoothetical 
graphite  epo.xy  laminate  approximating  5208-T-500  having  plies  1"3  microns 
thick  separated  by  9.90  micron  epoxy  boundary  layers.  The  velocity  in  the 
plys  is  taken  to  be  3300  m/s  and  the  velocity  in  the  epoxy  interface  2200  m, s . 
The  width  of  the  discontinuity  in  this  curve  is  roughly  .85  '■!Hz  as  determined 
from  tabulated  data,  while  the  approximation  (46)  gives  1.06  '■IHz  indicating 
that  an  error  of  about  20*!;  will  occur  when  this  formula  is  applied  to  measuring 
bond  line  thickness  in  graphite  epoxy.  The  discontinuity  in  C(jj)  may  also  be 
seen  in  a plot  of  w vs  '.p  as  seen  in  figure  5 for  various  values  of  ply  and 
bondline  thicknesses.  A graphical  measure  of  the  discontinuity  is  obtained  by 
extending  the  linear  portion  of  the  curve  into  the  non-travelling  wave  region 
and  measuring  the  frequency  displacement  between  the  two  lines.  As  is  shown 
by  table  II,  this  displacement  does  not  correlate  very  well  with  4f;  however 
it  does  provide  a number  which  is  easily  generated  from  empirical  data. 

EXPERIMENT 

Reflection  spectroscopy  measurements  (reflected  amplitude  vs  frequency) 
were  carried  out  on  a number  of  graphite  epoxy  panels  using  the  technique 
described  in  reference  1.  All  panels  were  prepared  and  cured  according  to 
manufacturers  specifications. 

RESULTS  .AND  DISCUSSION 
Experi.mental  Velocity  Measurements 

Figure  6 is  a typical  experimental  plot  of  Log^Q  of  amplitude  vs.  fre- 
quency for  ultransonic  waves  reflected  from  a graphite  epoxy  panel.  By 
measuring  the  positions  of  frequency  resonances  on  such  a plot  and  utilizing 
expression  (28)  it  is  possible  to  compute  dispersion  relationships  (velocity 
of  sound  vs.  frequency)  which  may  be  compared  with  analytical  predictions. 
Figure  ^ is  a plot  of  three  such  dispersion  relations  for  45°  graphite  epoxy 
laminates  of  4,  8 and  12  ply  thickness.  The  material  used  for  these  laminates 
was  5208-T-300  graphite  epoxy  with  a nominal  5 mil.  (Id"  micron)  ply  thickness. 

The  actual  values  of  the  laminate  thicknesses  are  shown  in  table  III  along 
with  the  resultant  average  ply  thicknesses.  These  values  were  used  in  the 
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1 "T"  ’ ~ 3 ^ 

V 

FIGURE  5.  Dispersion  Relations  For  Hypothetical  Laminates 
of  Varying  Thickness 


TKiUHL  6.  Ultrasonic  Keflectance  Spectrum  Of  21  Ply  .S208-'l-S00  (iraplute  lipoxy 


I'HJtlKli  7.  C()mi)ari son  OF  Lxj)er imeiita  1 And  llieoret  ical  KeFlectance 
Spectra  lor  (Iraph i t e/lipoxy  Laminates 
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tabu;  hi  - avhrace  ply  thicknesses  for  5208-T-300  graphite/epoxy  laminates 


LAMINATE 

TOTAL 

THICKNESS 

AVERAGE  PLY 
THICKNESS 

4 PLY  + 45° 

.021" 

.0052 

8 PLY  + 45° 

.043" 

.0054 

12  PLY  + 45° 

.062" 

.0052 

7 PLY  ♦ 0° 

.038" 

.0054 

For  the  materials  with  6 mil  nominal  ply,  thickness  variation  was  large  (>H0%) 
and  the  nominal  value  was  used  in  velocity  calculations. 
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velocity  calculations  for  5 mil  graphite  epoxy.  Ply  thicknesses  for  the  + 45° 
laminates  are  fairly  uniform  as  can  be  seen  in  the  micrograph  in  figure  8A. 

Frequencies  of  thickness  resonances  were  measured  to  about  ^ .04  Mllz  so 
that  corresponding  accuracies  of  velocity  measurements  ranged  from  ^ 4%  120 

m/s)  in  the  1 Mllz  region  to  ^ .4%  12  m/s)  near  10  Mllz. 

Comparison  of  Experimental  and  Theoretical  Dispersion  Relations 

Figures  7,  9 and  10  show  comparisons  of  theoretical  dispersion  relations 
(calculated  from  expressions  (34)  and  (29))  with  those  determined  from  ultra- 
sonic reflectance  measurements.  The  values  of  parameters  used  for  the  analyti- 
cal curves  and  their  sources  are  summarized  in  table  IV. 

Examination  of  these  experimental  and  theoretical  data  yields  the 
following  comparisons. 

1.  Locations  of  velocity  discontinuities  are  predicted  with 
reasonable  accuracy. 

2.  Experimentally  there  is  little  or  no  indication  of  the  predicted 
velocity  decrease  at  frequencies  just  below  a discontinuity. 

3.  The  functional  form  of  the  dispersion  relation  just  above  a 
discontinuity  is  described  reasonably  well  by  theory  but 
experimental  results  do  not  show  as  sharp  a discontinuity  as 
predicted . 

4.  Widths  of  discontinuities  are  rather  difficult  to  accurately 
measure  by  these  experiments  however  there  is  a definite 
indication  of  gaps  of  about  the  size  predicted. 

5.  A rapid  decrease  in  velocity  below  4 MHz  consistently  appears  in 
all  the  measurements  and  is  too  pronounced  to  be  attributed  to 
experimental  error.  It  is  probably  a result  of  material  disper- 
sion related  to  viscoelastic  damping. 

A number  of  reasons  can  be  hypothesized  for  the  lack  of  agreement  between 
theory  and  experiment,  the  least  important  of  which  is  probably  experimental 
error . 

It  is  clear,  from  figure  7 that  the  velocity  variation  among  samples  of 
the  same  material  can  be  significant  and  may  be  a function  of  laminate  thick- 
ness. Such  thickness  dependent  variations  would  of  course  contradict  results 
^ for  the  model  proposed. 

This  model,  in  fact,  has  a number  of  obvious  deficiencies: 

1 . Effects  of  material  dispersion  and  fiber  scattering  are  ignored 
and  these  effects  most  probably  contribute  strongly  to  the 
decrease  in  velocity  below  4MHz. 

2.  The  assumption  of  a perfectly  periodic  structure  can  often  be 
quite  bad  (this  can  be  seen  in  the  extreme  case  of  figure  8B) 
Variations  in  ply  thickness  would  undoubtedly  produce  the 
observed  lack  of  sharpness  in  the  velocity  discontinuity; 
however,  it  is  difficult  to  see  why  they  would  effect  the 
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FIGURE  8A.  Micrograph  Of  Cross  Section  Of  12  Ply  Graphite  Epox 
+ 45°  Laminate  (IlOx) 


FIGURE  8B.  Micrograph  Of  Cross  Section  Of  7 Ply  0°  Laminate 
Graphite  Epoxy  C112x) 
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Comparison  Of  experimental  And  Theoretical  Keflectance 
Spectra  lor  Crai)hite/epoxy  l.aminates 
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TABLE  IV  - PAR.AMETERS  USED  IN  THE  CALCULATION  OF  DISPERSION  RELATIONS 


PARAMETER 

VALUE 

SOURCE 

C.A 

3300  M/S 

t ULTRASONIC  GROUP  VELOCITY  ME.ASUREMENT 
FROM  5 MIL  LAMIN.ATE 

Cb 

2200  M/S 

COMPUTED  FROM  ESTIM.\TES  OF  ELASTIC  CONST.ANT 
.AND  DENSITY  C =\jCu/p 

da 

124u  (5  MIL) 
178u  (6  MIL) 

.975  .X  .NOMINAL  PLY  THICKNESS 
.9"5  X MA.X  THICKNESS  ME.ASUREMENT 

da/db 

40 

MICROSCOPY  ON  5 MIL  LAMINATE 

R 

. 535 

1 

ULTRASONIC  PULSE  REFLECTION 

This  is  believed  to  be  a better  estimate  of  the  velocity  of  sound  in  a layer 
than  can  be  calculated  from  static  values  of  modulus  and  density  since  the  latter 
calculation  takes  no  account  of  material  dispersion  and  gives  a value  of  about 
2300  M/S  which  is  clearly  to  low. 
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absence  of  a velocity  dip  which  was  predicted  just  below  the 
discontinuity. 

3.  Along  the  same  lines,  the  assumption  of  a periodic  laminate 

implies  the  existence  of  a resin  boundary  layer  of  well  defined 
thickness  on  the  surface  of  the  laminate.  Although  such  a layer 
undoubtedly  exists  its  thickness  probably  does  not  match  the  dg 
for  the  internal  boundary  layer.  Again  this  would  imply  that 
the  concept  of  a material  property,  namely  a thickness  indepen- 
dent phase  velocity  may  have  limited  validity  for  composites. 
This  topic  will  be  discussed  further  in  the  next  section. 

Variations  Among  Experimental  Results 

Measurements  of  sound  velocity,  in  addition  to  differing  from  theory  in 
functional  form  also  fail  to  exhibit  predicted  consistency  between  laminates 
of  different  thicknesses.  In  principal  curves  for  all  three  of  the  laminates 
in  figure  7 should  agree  to  within  experimental  error;  however  the  8 ply 
laminate  is  significantly  different  from  the  4 ply  and  12  ply  composites, 
inspite  of  the  fact  that  all  three  laminates  were  prepared  together  from  the 
same  starting  material. 

Even  slight  material  variation  can  cause  significant  changes  as  is  seen 
from  figure  9 which  is  the  dispersion  relation  for  the  same  material  as  shown 
in  figures  7 and  10  except  with  a nominal  6 mil  ply  thickness  instead  of  the 
standard  5 mils.  Clearly  the  functional  form  is  significantly  altered. 


SUMMARY 

Experimental  and  theoretical  efforts  applied  to  ultrasonic  NOT  of 
laminates  have  led  to  the  following  results; 

• A convenient  matrix  algebra  approach  can  be  constructed  for  treating 
problems  in  ultrasonic,  normal  incidence,  reflection  spectroscopy. 

• If  the  correct  boundary  conditions  are  imposed  a thickness  independent 
phase  velocity  can  be  defined  and  measured  for  periodic  laminates  in 

a manner  analogous  to  that  for  monolithic  materials. 

• Predicted  dispersion  relations^  for  finitely  thick  laminates  have 
discontinuities  which  are  the  same  as  those  for  infinite  laminates. 

• In  the  low  frequency  range  (<3MHz)  phase  velocity  is  constant  and 
dependent  only  on  laminate  composition  and  not  on  the  thickness  of 
specific  laminae. 

• Measured  effects  of  geometric  dispersion  in  graphite  epoxy  laminates 
are  in  semiquant itative  agreement  with  predictions.  Velocity 
discontinuities  are  on  the  order  of  1 MHz  both  theorically  and 
experimentally. 

• Measurable  variations  in  phase  velocity  are  greatest  at  high  frequencies 
(>8MIz)  and  are  less  that  about  6%. 

• Effects  of  geometric  dispersion  below  the  first  discontinuity  (-vSMHz) 
appear  weak  and  material  dispersion  dominates  resulting  in  variations 
in  phase  velocity  up  to  10%. 

• Because  of  variability  in  interface  and  lamina  thicknesses,  many 


-29- 


NADC-77093-30 


graphite  epoxy  laminates  are  only  a rough  approximation  of  the  idealized 
laminates  modeled  above.  As  a result  material  property  scatter  and 
thickness  dependent  material  properties  must  be  expected  in  real  material. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  use  of  spectrum  analysis  and  the  computation  of  dispersion  relations 
is  a very  sensitive  tool  for  interrogating  laminates,  in  that  effects  from 
interfaces  of  microns  in  thicknesses  can  be  revealed  with  acoustic  waves  of 
12MHz  and  below.  The  formalism  developed  here  may  be  useful  in  understanding 
such  spectra  in  appropriate  regimes. 

lor  example  the  velocity  of  sound  above  the  first  frequency  discontinuity 
could  be  used  as  a sensitive  monitor  of  laminate  uniformity  and  interface 
quality  providing  that  a sufficiently  high  level  of  material  uniformity  could 
be  achieved.  The  absence  of  this  uniformity  produces  scatter  in  material 
properties  and  is  the  greatest  obstacle  to  the  efficient  inspection  in 
composite  laminates. 

Laminates  of  potentially  high  uniformity  such  as  adhesively  bonded  sheet 
metal  laminates  are  excellent  candidates  for  inspection  by  the  above  techniques 
and  should  probably  be  considered  in  future  studies.  At  low  frequencies 
(i.e.,  wavelengths  much  longer  than  lamina  thickness)  static  moduli  can  be 
related  to  velocity  measurements  by  a simple  rule  of  mixtures  (40)  which  takes 
into  account  interface  phase  retardation.  This  technique  may  have  application 
to  directionally  solidified  eutectics. 

Because  material  dispersion  appears  to  have  such  a large  effect  at  low 
frequencies  (<3MHz)  mechanisms  for  this  dispersion  such  as  viscoelastic 
damping  should  be  studied. 

The  combined  theoretical  and  experimental  approach  described  here  relates 
to  discrete  velocity  measurements  and  therefore  has  limited  resolution.  Future 
investigations  may  require  the  use  of  continuous  wave  techniques  in  order  to 
permit  measurement  of  material  properties  continuously  as  a function  of 
frequency. 
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